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VERTICAL  BEAMFORMER  SIMULATION 


— *The  performance  of  an  all-digital  beamformer  for  a 
12-element  line  array  has  been  investigated  for  several  com- 
binations of  steering  angle,  frequency,  and  the  number  of  bits 
carried  at  various  points  in  the  beamforming  operation.  The 
results  are  compared  with  a completely  linear  system,  and 
with  a linear  system  with  discrete  delay  intervals  and  averaging 

between  taps._^ 

\ 

The  digital  system  is  shown  in  Fig.  1.  Each  hydro- 
phone feeds  a separate  4-bit  quantizer.  The  quantizer  charact- 
eristic is 

Q(x)  = Int<£  + 8)  , -7  < ^ * 7 (1) 

o a 

- 0 • I « -7 

- 15  . | * 7 

Where  a is  the  quantizer  level  spacing,  and  Int(x)  = largest 
integer  < x. 

Each  quantizer  feeds  a 4-bit  delay  line  with  taps 
spaced  t seconds  apart.  An  N-position  delay  tap  interpolator 
is  employed  to  provide  higher  resolution  for  the  delay  intervals. 
This  device  outputs  a linear  average  of  the  outputs  of  two  ad- 
jacent taps.  The  output  of  the  Nth  tap,  for  an  input  X(t),  is 

/x(t  - kT>  (N-n)  + X [(t-Oc+pTj-i  , 


providing  an  approximation  to  a true  delay  of  (k  + j|)  t seconds. 
Since  the  interpolation  process  multiplies  the  output  of  the 
quantizer  by  N and  thus  increases  the  number  of  bits  to 
4 + log2N,  we  divided  by  N'  to  truncate  the  output  to  the  re- 
quired value. 
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The  outputs  of  the  tap  interpolators  are  weighted  by 
suitable  integer  coefficients  and  summed  to  form  the  beam. 


Two  reference  systems  are  employed  to  investigate  the 
various  sources  of  degradation  in  the  digital  system.  One  is  a 
full  linear  system  in  which  the  quantizers  have  been  removed  and 
the  delay  lines  are  assumed  to  have  infinite  resolution  so  that 
tap  interpolation  is  unnecessary.  The  other  system  is  linear 
with  respect  to  amplitude,  but  the  delay  lines  have  a resolution 
of  t seconds  and  the  tap  interpolation  is  retained. 


In  order  to  steer  the  beam  to  an  angle  we  need  a 

delay  of 

(j  - 6.5)  | sin  p 

for  the  j element,  where  d is  the  distance  between  hydrophones, 
and  c is  the  velocity  of  sound.  The  tap  interpolator  provides 
an  approximation  to  a delay  of  [kj  + (nj/N)]T.  The  values  of 
kj  and  may  be  calculated  for  each  element  as  follows: 
for  one  half  of  the  array,  7 a:  j a:  12 

- I«[(j  - 6.5)  6 + .5] 

k)  ■ MtH  <3) 

”J  - h - 

The  kj  and  n^  for  the  other  half  of  the  array  are  found  by 
symmetry: 

for  1 s j * 6 


kj  “ -k13-j  _1 
nj  " N ' n13-j  , 


Letting  the  signal  be  S(t), 
signal  appearing  at  the  jth 


arriving  from  an  angle  a so  that  the 
hydrophone  is 
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Sj  [c  ” CJ  - 6.5)  | sin  a]  , 

we  nay  write  the  output  of  the  j**1  tap  interpolator  as 

f(N-n.)  q[s(t  - (j  * 6.5)  | sin  a + D.)1 

T.(t)  - Int{ 3 — 1 £ Li 

J 1 N' 

n.  Q[s(t  - (j  - 6.5)  § sin  a + D.  + t)1 


A similar  expression  describes  the  output  of  the  j 
tap  interpolator  in  an  otherwise  linear  system: 


Cj(t) 


(N-n^)  sj^t-  (j-6.5)~  sin  a+D^j  + n^sj^t-  (j-6.5)~  sino+Dj+tj 


The  output  from  the  jcn  element  of  the  full  linear  system  is 

Lj(t)  - S[t  - (j  - 6.5)  (§  sin  a - | sin  p)]  ^ . (6) 

1 N 

The  coefficients  in  Eq.  (5)  and  in  Eq.  (6)  are  introduced 

to  make  the  gain  of  the  linear  systems  the  same  as  that  of  the 
quantized  system. 

The  system  outputs  are 

12 


y(t)  - Int[< 


I (HJ  v*>)  +l 


for  the  quantized  system, 
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x(t)  - 


l vi  cj(t> 

1-1 


for  the  linear  system  with  tap  Interpolation;  and 


I wj  Lj(t) 


z(t) 


for  the  full  linear  system. 


The  Wj  are  Integer  shading  coefficients  and  M Is  the 
round-off  divisor  to  limit  the  output  of  the  quantized  system 
to  a specified  number  of  bits. 

There  are  three  sources  of  degradation  in  a system  of 

this  type: 

1.  Harmonic  distortion  due  to  quantization. 

2.  Distortion  of  the  beam  pattern  because  of  the 
finite  delay  resolution. 

3.  Frequency  distortion  from  the  tap  averaging. 

The  quantization  introduces  harmonic  distortion  which 
may  be  computed  as  follows: 

Let  p be  the  normalized  correlation  coefficient  be- 
yz 

tween  the  output  of  the  digital  system  and  the  output  of  the  full 

linear  system.  Let  a ^ and  a J"  be  the  variance  of  the  output  of 

y z 

the  quantized  and  linear  systems,  respectively.  Then  the  total 

signal  power  at  the  output  of  the  quantized  system  is  °y^Py8^ 

and  the  total  power  output,  is  qJ",  The  power  in  the  harmonics 
2 2 y 

is  o (l-p„_  ).  The  total  harmonic  distortion  is 

J 7* 
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% distortion  - 


harmonic  amplitude  __ 
signal  amplitude  x 100 


rl  - p 


x 100. 


We  will  not  consider  the  degradation  of  the  beam  pat- 
tern in  detail  here.  Instead,  we  shall  assume  that  four  point 
interpolation  and  a delay  of  60  u sec  are  adequate,  and  proceed 
to  investigate  the  amplitude  distortion  caused  by  tap  interpola- 
tion. The  output  of  a tap  interpolator  in  a linear  system  is 


where 


(N  - Hj)  sin  nit  + n^  sin  u)(t  + t)  - R sin  (uut  + <p) 


R -\lN2  - 4n (N-n)  sin2  ^ 


n 8 in  ist 


..  - -1  T n si: 

* tan  LIT-  n "4-  n 


COS  IDT. 


If  n - N/2,  then 


R - N cos  and 


_ _ OUT 

V ■ ■y*  * 


so  that  the  output  is 


N sin  £(t  + cos  . (13) 

The  delay  from  this  interpolation  is  exactly  half  way 
between  the  taps  on  the  delay  line  and  independent  of  frequency. 
The  amplitude  is  decreased  with  increasing  frequency  by  cos  njr. 
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For  N » 4,  n ■ 1 or  3 we  find  for  the  amplitude  function 


R - 2 \ 4 - 3 sin2  4jf 


For  the  delay,  with  n - 1 


cp  ■ tan"  I' 


-If  sin  hit 


cos  iutJ 


and  for  n - 3 


cp  ■ tan" 


-If  3 sin  w t y 
11  + 3 cos  uitj 


For  t ■ 60  usee  delays  are  shown  in  Table  I,  and  the  amplitudes, 
in  dB  relative  to  the  ideal,  are  shown  in  Table  II. 


TABLE  I 


TABLE  II 


2 kHz  3 kHz  Ideal 


2 kHz  3 kHz 


-0.47  -1.05 


-0.64  -1.48 


-1.05 
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l The  variations  in  effective  time  delay  and  shading 

t 

I with  frequency  will  introduce  variations  of  the  beam  pattern  with 

frequency  in  addition  to  those  caused  by  the  changing  ratio  of 
| wavelength  to  array  dimensions. 

\ The  equations  for  the  output  of  the  various  systems 

have  been  programmed  for  evaluation  by  a computer , using  a sine 
wave  signal.  The  program  calculated  the  correlation  coefficients 
and  relative  gains  of  the  digital  system  vs.  the  full  linear 
system  and  vs.  the  linear  system  with  tap  interpolation. 

The  results  are  summarized  in  Table  III,  and  plots  of 
the  output  of  the  digital  system  vs.  the  full  linear  system  are 
shown  in  Figs.  2 through  14. 

DISCUSSION  OF  RESULTS 

The  output  distortion  depends  on  the  steering  angle, 
the  frequency,  the  number  of  bits  in  the  tap  interpolator,  and 
the  shading  function.  If  we  use  total  harmonic  distortion  as  a 
criterion,  the  worst  cases  with  respect  to  steering  angle  are 
those  which  require  the  same  tap  on  the  tap  interpolators  for 
[ each  element.  For  instance,  at  0°  (perpendicular  to  the  column) 

the  output  of  each  element  is  taken  directly  from  a delay  line 
[ tap.  At  26.8°,  each  output  is  from  the  center  tap  of  the  tap 

interpolator.  These  are  worst  cases  since  all  of  the  quantizers 
step  at  the  same  time,  producing  a transfer  function  identical 
to  that  of  the  basic  quantizer.  For  other  steering  angles,  the 
output  is  the  sum  of  signals  taken  from  various  taps  on  the  inter- 
polators. This,  together  with  the  residual  phase  errors,  tends 
to  increase  the  effective  amplitude  resolution  of  the  system  and 
reduce  the  distortion.  The  relative  phase  error  increases  with 
frequency  so  that,  within  limits,  the  distortion  is  reduced  with 
increasing  frequency. 
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When  the  beam  is  steered  to  26.8°  from  broadside,  the 
output  of  each  element  is  taken  from  the  middle  tap  of  the  tap 
interpolator.  The  transfer  function  becomes  that  of  one  quantizer 
followed  by  the  tap  interpolator.  This  transfer  function  does 
not  have  a step  at  the  origin.  This  may  lead  to  difficulties 
with  weak  signal  suppression  and  needs  to  be  investigated  further. 

The  four  point  tap  interpolator  operating  on  the  out- 
put of  a four  bit  quantizer  requires  six  bits  to  maintain  the 
full  accuracy  of  the  system.  Truncating  the  output  to  four  bits 
increases  the  distortion.  However,  only  in  the  case  of  steering 
to  26.8®  does  the  distortion  rise  above  that  obtained  at  broad- 
side. Truncation  does  destroy  the  symmetry  of  the  transfer 
function.  The  quantizer  itself  will  produce  only  odd  harmonics. 
The  truncated  output  of  the  tap  interpolator  will  contain  even 
harmonics  as  well. 

It  should  be  noted  that  the  distortion  is  also  a 
function  of  the  sum  of  the  weights  used  for  shading.  This  number 
affects  the  way  in  which  the  final  output  is  resolved  after 
round-off.  For  the  particular  case  considered  here,  the  sum  of 
the  weight  is  72.  The  amplitude  of  one  output  step  at  broadside 
before  round-off  is  72  x 4 ■ 288.  The  maximum  output  is 
72  x 4 x 15  - 4320.  If  we  truncate  this  to  approximately  7 bits 
by  dividing  by  32,  we  find  that  the  maximum  output  is  133  and 
each  step  is  9.  This  is  resolved  exactly  at  the  final  output. 
However,  the  maximum  output  is  slightly  greater  than  can  be 
handled  by  7 bits. 

If  we  reduce  the  sum  of  the  weights  to  68,  each  step 
will  be  68  x 4/32  ■ 8.5  and  the  maximum  output  is  127.5.  Thus, 
we  find  that  the  broadside  output  has  alternate  steps  of  8 and  9 
units,  but  the  output  is  within  the  range  of  7 bits.  The  distor- 
tion is  increased  from  6.00%  to  6.077o  by  the  change  in  shading 
pattern. 
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Transfer  Function  for  Beam  Steered  to  8°  at  3 kHz  and  6 bit 
Tap  Interpolation 


- Transfer  Function  for  B< 
and  4 bit  Tap  Interpola 
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Fig.  6 - Transfer  Function  for 
2 kHz  and  4 bit  Tap  Ii 


Transfer  Function  for  Beam  Steered  to  20°  at 
3 kHz  and  6 bit  Tap  Interpolation 
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Transfer  Function  for  Beam  Steered  to  20  at 
3 kHz  and  4 bit  Tap  Interpolation 
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